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ABSTRACT 



Diabetes mellitus affects the brain. Both type 1 and type 2 diabetic patients are associated with white matter 
(WM) damage observable to diffusion tensor imaging (DTI). The underlying histopathological mechanisms, 
however, are poorly understood. The objectives of this study are 1) to determine whether streptozotocin 
(STZ)-induced type 1 diabetes is associated with WM damage observable to DTI; and 2) to understand the 
pathophysiological aspects underlying STZ-induced brain injuries. Male Sprague-Dawley rats received a single in- 
traperitoneal injection of STZ (62 mg/kg). DTI was used to assess brain abnormalities at 4 weeks after induction, 
combined with histological assessments and ultrastructural analysis. Compared to controls, the STZ-induced rats 
showed significantly reduced fractional anisotropy (FA) in the motor/somatosensory cortex and striatum. Histo- 
logically, the cortex and striatum of the diabetic animals are characterized by demyelination and axonal degrada- 
tion. In conclusion, STZ-induced diabetes is associated with striatal/cortical injuries observable to DTI. The DTI 
abnormalities are likely manifestations of demyelination and axonal degradation in the affected brain regions, 
and can potentially be used as surrogates for evaluating diabetic brain injuries. 

© 2012 The Authors. Published by Elsevier Inc. All rights reserved. 
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Introduction 

Patients with diabetes mellitus, and in particular those with poor 
glycemic control, often experience brain complications (van Harten 
et al., 2006). It is well known that long-standing diabetes increases 
the risk of brain atrophy, lacunar infarcts and white matter (WM) le- 
sions. The functional/behavioral consequences of diabetic brain com- 
plications include cognitive dysfunction and movement disorders 
(Shan et al., 1998; van Elderen et al, 2010; van Harten et al, 2006). 

Diffusion tensor imaging (DTI) is a magnetic resonance imaging 
(MRI) technique that has been increasingly used in the diagnosis and 
prognosis of WM diseases (Tournier et al., 2011). Parametric diffusion 
indices derived from DTI, such as fractional anisotropy (FA) and direc- 
tional diffusivity are often used to assess microstructural changes of 
brain tissues (Tournier et al., 201 1 ). 

Only few previous studies have assessed WM damage in diabetic 
patients with DTI. It was shown that patients with long-standing 



* This is an open-access article distributed under the terms of the Creative Commons 
Attribution-NonCommercial-ShareAlike License, which permits non-commercial use, 
distribution, and reproduction in any medium, provided the original author and source 
are credited. 

* Corresponding author at: Wuhan Institute of Physics and Mathematics, Chinese 
Academy of Sciences, 30# Xiaohongshan, Wuhan, Hubei 430071 , China. Tel.: + 86 27 
87198542; fax: +86 27 87199291. 

E-mail address: leihao@wipm.ac.cn (H. Lei). 
1 Mingming Huang and Lifeng Gao contributed equally. 



type 1 diabetes had significantly reduced FA in the posterior corona 
radiata and optic radiation, which correlated with the duration of di- 
abetes (Kodl et al., 2008). These patients also had significantly re- 
duced cortical thickness in brain regions with high connectivity to 
the injured WM tracts (Franc et al., 2011). DTI assessment among in- 
dividuals with type 2 diabetes demonstrated diffuse but predomi- 
nantly frontal and temporal WM microstructural abnormalities (Yau 
et al., 2009). A more recent DTI study found that, compared to the 
control, type 2 diabetic patients had significantly increased mean diffu- 
sivity (MD) of the brain parenchyma and significantly decreased FA in 
the bilateral frontal WM (Hsu et al., 2012). They further showed that 
there was a significant association between disease duration and micro- 
structural properties in the bilateral cerebellum, temporal lobe WM, 
right caudate and bilateral cingulate gyrus (Hsu et al., 2012). 

In a mouse model of streptozotocin (STZ) -induced type 1 diabetes, 
WM changes and brain atrophy were detected 8 months after induction 
(Francis et al., 2008). Yang et al. showed that T 2 -weighted and magne- 
tization transfer MR1 can be used to assess WM abnormalities in rats 
with type 1 diabetes or type 2 diabetes, and the changes of MR1 signal 
intensity are associated with loss of myelin proteins and oligodendroglia 
(Yanget al., 2011). 

In this study, we used DTI to assess brain abnormalities in a rat model 
of STZ-induced type 1 diabetes at 4 weeks after induction. The imaging 
results were correlated to those of histological assessments and ultra- 
structural analysis. The main purposes of the study are 1 ) to determine, 
without glycemic control, whether STZ-induced type 1 diabetes can 
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induce WM damage observable to DTI; and 2) to understand the patho- 
physiological aspects underlying STZ-induced brain injuries. 

Materials and methods 

Animals 

All animal protocols were approved by the Institutional Animal Care 
and Use Committee. Forty male Sprague-Dawley rats, about 10 weeks 
old and weighing about 300 g, were obtained from the Experimental 
Animal Center, Wuhan University. Type 1 diabetes (n = 24) was induced 
by a single intraperitoneal injection of streptozotocin (62 mg/kg) (Sigma, 
St. Louis, MO) dissolved in 0.01 mol/L citric acid solution (pH = 4.5). 
The control animals (n = 16) were injected with the same amount of 
solvent. The animals were then housed individually in a temperature- 
and humidity-controlled environment on a 12-h light/dark cycle. To 
prevent hypoglycemia, the STZ-induced rats were allowed to have 
free access to D5W (water with 5% dextrose) during the first 24 h 
after STZ injection. Three days after the induction, fasting blood glucose 
levels and body weights were measured. Six STZ-induced rats that had a 
fasting blood glucose level less than 13.5 mmol/L were excluded. Fasting 
blood glucose level and body weight were measured weekly for the next 
4 weeks. 

MRI acquisition 

All magnetic resonance experiments were performed on a 7 T/ 
20-cm-diameter bore Bruker Biospec scanner at 4 weeks after STZ in- 
duction. A volume coil was used for radiofrequency transmission, and a 
quadrature surface coil for signal detection. DTI was performed with a 
4-shot spin-echo echo planar imaging sequence, an encoding scheme 
of 30 gradient directions homogenously distributed on the unit sphere 
and the following parameters: repetition time (TR) 5000 ms, echo time 
(TE) 26 ms, field of view (FOV) 3 cm x 3 cm, slice thickness 0.8 mm, ma- 
trix size 128 x 128, A 14 ms, 6 3 ms, two b values of 0 and 800 s/mm 2 
and 4 averages. The animals were anesthetized by 1.5-2% isoflurane in 
pure 0 2 delivered via a nose cone during the imaging sessions. 

Histology and immunohistochemistry 

After being imaged, three randomly selected rats from each group 
were anesthetized with an overdose of 5% chloral hydrate, and per- 
fused transcardially with 0.9% saline and 4% paraformaldehyde in 
phosphate-buffered solution (PBS). The brains were immediately re- 
moved and postfixed in the same fixative at 4 °C until being sectioned 
on a ciyostat (Leica, Germany). Coronal brain sections 10 um-thick were 
obtained and stored at — 20 °C until use. 

Bielschowsky's sliver (BS) staining was used to assess the morpholo- 
gy and integrity of axons and neuropil. The brain sections were washed 
with distilled water and incubated in 10% silver nitrate solution for 
35-40 min at room temperature. The sections were then processed se- 
quentially with a 5% methanal solution for 5 min, a 1 0% ethanol solution 
for 2 min, a silver amino solution for 2 min and a 10% methanal solution 
for 5 min, followed by washes with distilled water and dehydration in 
100% ethanol. 

Kliiver-Barrera (KB) staining was used for demonstrating the exis- 
tence of myelin. The brain sections were incubated with a 0.1% Luxol 
fast blue (LFB) solution at 60 °C overnight and rinsed with a 95% etha- 
nol solution and distilled water. The brain sections were then immersed 
in a 0.05% lithium carbonate solution for 1 min and a 70% ethanol solu- 
tion for 2 min, and treated with 0.1% cresyl violet and rinsed with 95% 
ethanol. Finally, the sections were dehydrated in 100% ethanol, cleared 
with xylene and mounted. 

Nissl staining was used to assess neuronal damage in the cortex. The 
brain sections were incubated with a 5% toluidine blue solution at room 



temperature for 1 5 min. Following rinses with tap water, the sections 
were dehydrated and mounted. 

Diminution in neurofilament immunostaining of SMI-31 is usually 
considered a sign of axonal degeneration (Budde et al., 2007). For 
SMI-31 immunohistochemistry, the brain sections were pre-incubated 
with a 0.01 M citrate buffer (pH 6.0) at 80-85 °C for about 10 min, 
until the temperature of the antigen retrieval solution had reached 
room temperature. After a brief wash in 0.01 M PBS, the brain sections 
were incubated in 0.2% Triton X-100 at 37 °C for 30 min, washed in 
0.01 M PBS and then incubated with 10% normal goat serum at 37 °C 
for 40 min. Subsequently, the brain sections were incubated in anti- 
phosphorylated neurofilaments SMI-31 antibody (diluted 1:1000 in 
PBS; Stemberger Monoclonals Lutherville, Maryland) at 37 °C for 1 h 
and then at 4 °C for additional 21 h. Further washes with 0.01 M PBS 
were performed before incubation in Alexa 488-rabbit anti-goat IgG an- 
tibody (diluted 1 :200 in 0.01 M PBS) at 37 °C for 1 h. The sections were 
washed with 0.01 M PBS and dehydrated before being mounted. 

Methylene blue staining was used for morphometric analysis of 
myelinated axons. Transmission electron microscopy was used for ultra- 
structural assessments of myelin integrity and axonal injuries. Samples 
from the striatum and cortex were collected and fixed with 2.5 wt.% 
glutaraldehyde dissolved in 0.01 M PBS (pH 7.4). The samples were 
postfixed with 1% osmium tetroxide, dehydrated, embedded in resin 
and sectioned to 80 nm- or 1 pm-thick slices. The 1 urn-thick slices 
were stained with methylene blue. The 80 nm-thick slices were stained 
with uranyl acetate and lead citrate for ultrastructural tissue examination 
under an H-7000FA transmission electron microscope (Hitachi, Japan). 

Image pre-processing 

The DTI data from two rats in the STZ-induced group were excluded 
because of poor image quality. As a result, the DTI data from a total of 
sixteen rats in the STZ group and sixteen rats in the control group 
were included in the final analysis. For each animal, the FMRIB's Diffu- 
sion Toolbox in FSL (http://www.fmrib.ox.ac.uk/fsl) was used to realign 
diffusion-weighted images to corresponding non-diffusion-weighted 
(b = 0) images with an affine transformation to minimize image distor- 
tion due to eddy currents. The diffusion tensor for each voxel was esti- 
mated by a multivariate linear fitting algorithm, and then diagonalized 
to obtain its eigenvalues (\ lt \ 2 , \a). The voxel-wise values of FA, axial 
diffusivity D a ) and radial diffusivity ( (\ 2 + X.3)/2, D r ) were calculat- 
ed. Residual maps were calculated as previously described to evaluate 
the accuracy of tensor fitting (Tournier et al., 201 1 ). 

Voxel-based analysis of DTI data 

Voxel-based analysis (VBA) of FA data was performed using statistical 
parametric mapping (SPM8, http://www.fil.ion.ucl.ac.uk/spm/). The fol- 
lowing processing schemes were used: (1) creating population-specific 
template. First, the b = 0 image of one representative control rat was 
trilinearly interpolated into an isotropic voxel size of 0.117 mm 3 and 
used as the initial template. Then the b = 0 images of all the other 
(including control and STZ-induced) animals were co-registered to 
this initial template using affine transformation. A population- 
specific template was then created by averaging the normalized images; 
(2) VBA on FA. The b = 0 image of each rat was co-registered to the 
population-specific template using affine transformation. For each indi- 
vidual animal, the same transformation matrix was applied to normalize 
the FA maps spatially. All normalized FA images were then smoothed 
with an isotropic 0.351 -mm (3 times the dimension of the isotropic 
voxel after interpolation) full width at half maximum (FWHM) Gaussian 
kernel to reduce the effect of misregistration in spatial normalization. 
The normalized and smoothed FA images were then analyzed voxel- 
wise by SPM8 within the framework of general linear model. Only the 
pixels where all rats had a FA value greater than 0.15 were analyzed. 
This empirical cut-off value was used to eliminate fictitious results arising 



M. Huang et al I Neurolmage: Clinical I (2012) 57-65 



59 



from inter-subject variations and in the voxels with a low probability of 
being WM, such as the CSF voxels (Hugenschmidt et al., 2008). The t 
maps were thresholded using a confidence level ofp<0.005 (t>3) and 
an extent cluster size greater than 10 contiguous voxels. 

Histogram and region of interest (ROI) analysis of DTI data 

The DTI changes in the striatum and cortex were further analyzed 
by histogram and ROI analysis as previously described (Yu et al., 
2006). Masks of the bilateral striatum and cortex (Figs. 2 and 3) 
were derived manually from the FA template. The ROI representing 
bilateral striatum spanned 24 interpolated slices from Bregma 2.52 
to 0. The ROI representing bilateral cortex spanned 43 interpolated 
slices from Bregma 3.00 to —0.60. The resulting masks were then 
superimposed onto the FA images of each individual animal. Histo- 
grams containing 1000 bins were created from the bilateral striatum 
and cortex of each rat. To compensate for variability of brain size, each 
bin was normalized by the total number of voxels contributing to the 
histogram. To reduce the possible effects of local minimum, the histo- 
grams from each rat were subjected to a 5-point median smoothing 
procedure. Three parameters were then derived from the smoothed 
histogram of each subject, including the averaged FA value, histogram 
peak height and histogram peak location. Group-averaged histograms 
were also calculated for presentation. 

The striatum was composed of fiber-rich WM-like regions having 
high FA values and gray matter (GM)-like regions having low FA values. 
To pinpoint the exact locations of the diabetes-induced FA changes, we 
divided the ROI representing the whole striatum into two sub-ROIs 
based on the averaged FA map from the control group. A cutting FA 
value of 0.2 was used to classify each voxel within the striatal mask 
into either the high FA sub-ROI or the low FA sub-ROl (Smith et al., 
2006). The averaged DTI indices within the sub-ROIs were calculated 
and compared between the control and STZ groups (Fig. 3). Similar 
analyses were also conducted for the cortex (Fig. 4). 

Statistical analysis 

Repeated measures analysis of variance (ANOVA) was used to ana- 
lyze body weight and blood glucose data. Two-tailed Student's t-tests 
were used to assess the statistical significance of inter-group difference 
in histogram parameters and diffusion indices. The significance level 
was set at p< 0.05 with Bonferroni correction for multiple comparisons. 

Results 

Body weight and blood glucose level 

Table 1 shows the body weight and fasting blood glucose data. No 
statistically significant difference in body weight was found between 
the groups before saline/STZ treatment, nor in fasting blood glucose 
level. The control rats showed increasing body weight over time, 
while the STZ group had decreasing body weight. The fasting blood 
glucose level in the STZ group was significantly higher than that in the 
control group at 3 days and the time points after. Repeated measures 
ANOVA revealed that the main effects of group, time and group x time 



interaction were statistically significant for both body weight and fasting 
blood glucose level. 

Voxel-based analysis of DTI data 

Fig. 1A shows representative raw diffusion-weighted (b = 0 and 
800 s/mm 2 ) images acquired from a control animal, as well as a residual 
map demonstrating the accuracy of tensor fitting. The results of VBA are 
shown in Fig. IB. Compared to the control animals, the diabetic rats 
exhibited significant decreases of FA in the bilateral striatum, motor cor- 
tex and somatosensory cortex. 

Histogram/region of interest analysis of DTI data 

Compared to the control group, the STZ group showed a left-shifting 
FA histogram in the bilateral striatum (Fig. 2A). The inter-group differ- 
ence in histogram peak location (Fig. 2B) was statistically significant, 
but not the difference in histogram peak height (Fig. 2C). The STZ 
group showed significantly decreased average striatal FA, which oc- 
curred in both the GM-like regions where FA < 0.2 and the fiber 
bundle-dominated region where FA>0.2 (Fig. 2D). The reduction 
of striatal FA in the STZ group was accompanied by a significantly de- 
creased D a and an unchanged D r (Fig. 2E and F). 

The group-average FA histograms of the bilateral cortex are shown 
in Fig. 3A. Compared to the control group, the STZ group had a signifi- 
cantly (p<0.05) higher histogram peak height (Fig. 3B), but unchanged 
histogram peak location (Fig. 3C). The decrease of cortical FA in the STZ 
group appeared to be confined to the WM-like regions where FA>0.2 
(including the molecular layer and the pyramidal layers, Fig. 3D and E). 
The changes of FA in the regions with FA < 0.2 (including the granular 
layers and part of the polymorphic layer) were not statistically significant 
(Fig. 3D). For the STZ group, the decrease of FA in the cortical region with 
FA > 0.2 was accompanied by a slightly increased D r and an unchanged D a 
(Fig.3F). 

Histopathological changes in striatum 

The histopathological changes of the fiber bundles in the striatum 
of STZ-induced rats included: less staining with LFB (Fig. 4B) and 
SM1-31 (Fig. 4D), and sparser/disarranged neurofibrils as revealed 
by BS staining (Fig. 4F). The result of both BS staining and methylene 
blue staining (arrows in Fig. 4H) indicated vacuole formation in the 
striatal fiber bundles of the STZ-induced animals. Ultrastructural 
analysis of the axons in the striatal fiber bundles showed signs of rar- 
efaction of myelin sheath (arrows in Fig. 4J), myelin loss (arrowhead 
in Fig. 4K) and axonal degradation (asterisk in Fig. 4I<). 

Histopathological changes in cortex 

The histopathological changes in cortex of the STZ-induced ani- 
mals included rarefied neuropil in the molecular layer (Fig. 5C), dys- 
trophic pyramidal neurons exhibiting fewer processes and poorly 
BS-stained somas (Fig. 5E and G), and signs of neuronal dystrophia 
(arrows in Fig. 50) or loss (asterisk in Fig. 5M) in the polymorphic 
layer. Ultrastructural analysis showed that the pyramidal neurons of 
the STZ animals were associated with compromised integrity of the 



Table 1 

Fasting blood glucose level and body weight in control (Con) and STZ-treated animals. 







DayO 


Day 3 


Day 7 


Day 14 


Day 28 


Blood glucose (mmol/L) 
Body weight (g) 


Con (n = 16) 
STZ (n = 16) 
Con (n = 16) 
STZ (n = 16) 


4.2 ±0.6 
4.0 ±0.5 
293.3 ±27.5 
300.7 ±17.3 


5.4 ±1.1 
23.6 ±7.6* 


3.5 ±0.5 
19.0 ±6.3* 
322.2±31.0 
268.0±21.7* 


3.4 ±0.3 
17.3 ±3.9* 
343.5 ±30.4 
250.8±31.5* 


3.7 ±0.6 
25.9 ±7.0* 
387.6 ±36.1 
230.4 ±30.2* 



* Significantly different from the control group (p<0.0001, two-sample t-test). 
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Fig. 1. Representative raw diffusion-weighted images from a control rat are shown in (A), as well as the corresponding residual map demonstrating the accuracy of tensor fitting. 
Panel B shows that the STZ-induced rats had diffuse reduction of FA in bilateral striatum and motor/somatosensory cortex. Pixels with statistical significance (t>3. voxel extent 
threshold = 10) are shown in color and superimposed on the population-specific FA template generated from 16 control rats and 16 STZ-induced rats. An FA threshold of 0.15 was 
used when generating the t map. The color bar in (A) denotes the relative amplitude of the residuals (as percentage of the signal intensity of the corresponding pixel in the raw images). 
The color bar in (B) denotes t scores (range 0-8). 




Whole FA > 0.2 FA < 0.2 Da Dr Da Dr 



Fig. 2. Histogram and region of interest analysis of the diffusion indices in the bilateral striatum. The group-averaged FA histograms of the bilateral striatum (from Bregma 2.52 to 
0 mm) are shown in the form of line graphs in A. Compared to the control group (Con), the STZ group (STZ) showed a statistically significant (*p<0.001 ) reduction in histogram 
peak location (B), but no change in histogram peak height (C). Compared to the control group, the STZ group showed significantly reduced average FA in the whole striatum, the 
gray matter-like region with FA<0.2, and the white matter-like region with FA>0.2 (D, *p<0.001 ). The reduction of FA in the diabetic striatum was associated with a significantly 
decreased axial diffusivity (D a , F and G, *p<0.05 vs. Con) and unchanged radial diffusivity (D r ). 
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Fig. 3. Histogram and region of interest analysis of the diffusion indices in the bilateral cortex The group-averaged FA histograms of the bilateral cortex (from Bregma 3.00 to — 0.60 mm) 
are shown in the form of line graphs in A. Compared to the control group (Con), the STZ group (STZ) showed a statistically significant (*p<0.05) increase in histogram peak height (B), but 
no change in histogram peak location (C). The STZ group showed significantly reduced FAonly in the cortical region with FA> 0.2 (D and E). Panel D shows the regions with FA> 0.2 (in red 
color) within the cortex mask. The FA changes in the whole cortex and cortical region with FA< 0.2 (including granular layer and partly polymorphic layer) were not statistically significant 
(E). Panel Fdemonstrates that the reduction of FA in the cortical region with FA> 0.2 was associated with a slightly increased radial diffusivity (D r ) and an unchanged axial diffusivity (D a ). 



myelin sheath surrounding the axons (arrows in Fig. 5I), loss and/or 
fragmentation of neurofilaments in the dendrites (asterisk in Fig. 5I<), 
mitochondrial abnormalities (arrowheads in Fig. 51 and K) and aggrega- 
tion of cellular organelle (circle in Fig. 5K). 

Discussion 

In the present study, brain injuries in a rat model of STZ-induced 
type 1 diabetes were evaluated with DTI. It was shown that, without 
glycemic control, the STZ-induced animals developed reduced FA in 
the striatum and motor/somatosensory cortex as early as 4 weeks 
after induction. Histologically, the brain regions showing FA changes 
were associated with demyelination, axonal degradation and neuro- 
nal dysfunction/loss. 

Diabetes-induced DTI changes 

Diffusion indices derived from DTI, such as FA, are widely-used 
surrogates for WM microstructural integrity. Clinical studies have 
shown that both type 1 diabetic and type 2 diabetic patients are associ- 
ated with WM abnormalities observable to DTI. Lower FA values were 
found in the posterior corona radiata and the optic radiation of type 1 
diabetic patients (Kodl et al., 2008), and in the frontal, temporal and 
right caudate WM of type 2 diabetic patients (Hsu et al., 2012; Yau 
et al, 2009). 

Interestingly, we did not find any significant FA changes in the major 
WM tracts of diabetic animals. Instead, without glycemic control, the 
STZ-induced rats were shown to have significantly reduced FA in the bi- 
lateral striatum and motor/somatosensory cortex as early as 4 weeks 
after induction. The FA reduction in the striatum appeared to be homog- 
enous, and accompanied by a significantly decreased D a and a slightly in- 
creased D r (Fig. 2). It is generally believed that the increase of D r in WM is 
associated with demyelination (Song et al., 2002, 2005). In contrast, D a 



appeared to be a marker more closely related to axonal damage 
(Sun et al., 2006). Decreased D a has been shown to be associated with 
axonal degeneration (i.e., loss of phosphorylated neurofilaments) and 
impaired axonal transport (Deboy et al., 2007; Wu et al., 2007). Histo- 
logically, the striatal fiber bundles in our STZ-induced animals were 
characterized by both demyelination and axonal degeneration, consis- 
tent with the DTI findings. 

The water diffusion anisotropy in the cortex mainly originates from 
the radial pyramidal neurons and the packed arrangement of axons 
(Mukherjee et al., 2002). It has also been reported that neuron density 
and glial activity can also modulate anisotropic diffusion of water in the 
cortex (Beppu et al., 2005; Neil et al., 1998). Unlike the case in the stria- 
tum, the cortical changes of FA were not homogenous, but rather con- 
fined to layers with high FA values (including the molecular layer, the 
pyramidal layer and part of the polymorphic layer, Figs. IB and 3). The re- 
duction of FA in the cortex was mainly caused by an increase in D r , with- 
out apparent changes in D a (Fig. 3 ). Interestingly, the cortical regions with 
FA reduction appeared to coincide with those showing axonal degrada- 
tion and/or neuronal loss. For example, the neuropil in the molecular 
layer of the STZ-induced rat appeared to be rarefied compared to that 
in the control; the neurons in the pyramidal layer of the STZ-induced an- 
imals exhibited compromised integrity of the myelin sheath, damaged 
neurofilaments in the dendrites and mitochondrial abnormality (Fig. 5). 
Therefore degeneration of the pyramidal neurons could have played a 
significant role in causing the diffusion changes in the cortex. 

However, it should be emphasized that the rat striatum and cortex 
are heterogenous regions in which WM and CM are mixed together. 
Due to the complex tissue configurations in these regions and the 
partial volume effects, it is the most likely that there is more than 
one dominant fiber orientation within a single striatal/cortical voxel. 
Under such circumstances, the calculation of parametric diffusion in- 
dices may not be as straightforward as that in major WM tracts, and 
could be affected by factors that are not necessarily related to WM 
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Fig. 4. Histopathological staining of the striatal sections from control rats (Con, panels A, C, E, Gand I) and STZ-induced rats (STZ, panels B, D, F, H, J and K), including KB staining (Aand B), 
SMI immunohistochemistry (C and D), BS staining (E and F), methylene blue staining (G and H) and uranyl acetate-lead citrate staining (I to K). Compared to that of the control rat, the 
striatal fiber bundles in the STZ-induced rats are less intensively stained with LFB (B) and SM1-31 (D), and associated with disarranged neurofibrils and vacuole formation (arrows in Fand H). 
Ultrastructural analysis showed that the striatal fiber bundles of the STZ-induced animals are characterized by rarefied myelin sheath (arrows in J), myelin loss (arrow head in K) and signs of 
axonal degradation (asterisk in K). 



integrity (Jeurissen et al., in press; Vos et al., 2012; Wheeler-Kingshott 
and Cercignani, 2009). Caution thus should be exercised in linking 
diabetic-induced FA, D a and D r changes to specific histopathological 
mechanisms. 

In contrast to the observations in the rat model of diabetic ketoacidosis 
(Lam et al., 2005), we found no significant reduction of cortical apparent 
diffusion coefficient (ADC) in the diabetic rats (data not shown). Cere- 
bral edema is the most common complication of diabetic ketoacidosis 
(Glaser et al., 2001), which may have contributed to the ADC increase 
observed by Lam et al. (2005). In this study, no evident signs of cerebral 
edema were observed in the diabetic animals. 

Histological changes in diabetic striatum and cortex 

The histological changes observed in the striatum and cortex of the 
diabetic animals are consistent with the findings from previous studies, 



showing that the diabetic brain is associated with demyelination and 
axonal injury. For example, smaller cells with reduced myelin content 
were found in the cortex of the STZ-induced diabetic animals at 
4 weeks (Malone et al., 2006). Hernandez-Fonseca et al. showed 
that disarrangement of myelin sheath and oligodendrocyte abnor- 
malities can be observed in STZ-induced rats at 8 weeks after induc- 
tion (Hernandez-Fonseca et al., 2009). Reductions in myelin quantity 
and loss of oligodendrocytes were observed in both type 1 diabetic rats 
and type 2 diabetic rats at 4 months (Yang et al., 2011). In a mice 
model of STZ-induced type 1 diabetes, myelin loss occurred in both 
WM and CM regions at 8 months after induction (Francis et al., 2008). 
On the other hand, impairment of axonal transport in the olfactory 
bulbs of the STZ-induced mice occurs as early as one week after induction 
(Sharma et al., 2010). It was demonstrated that, after 8 weeks of hyper- 
glycemia, the pyramidal cells in the parietal cortex, medial prefrontal cor- 
tex and hippocampus of the STZ-induced diabetic rats showed dendritic 



M. Huang et al I Neurolmage: Clinical I (2012) 57-65 63 




Fig. 5. Histopathological staining of the cortical sections from control rats (Con, panels B, D, F, H, J, Land N) and STZ-induced rats (STZ, panels A, C, E, G, I, K, M and 0), including BS 
staining (A to E), methylene blue staining (F, G, N and 0), uranyl acetate-lead citrate staining (H to K) and Nissl staining (L and M). Panel A illustrates the laminar structures of the 
cortex. The neuropil in the molecular layer of the STZ-induced rat (C) appeared to be rarefied compared to that in control (B). The arrows in D, F and N indicate normal-appearing 
neurons in the pyramidal layer and polymorphic layer. The arrow in E indicates the poorly BS-stained soma of a pyramidal neuron, while those in G and 0 indicate dystrophic neu- 
rons. The results of ultrastructural analysis were presented in H-K, in which arrowheads indicate mitochondria, the letter d indicates dendrites, and asterisks indicate 
neurofilaments in the dendrites. Compared to the control animals, the neurons in the pyramidal layer of the STZ-induced animals exhibited compromised integrity of the myelin 
sheath (arrows in I), mitochondrial abnormality (arrowheads inland K), loss and/or fragmentation of neurofilaments in the dendrites (K), and cellular organelle aggregation (circle in K). 
The asterisks in M indicate regions with evident neuronal loss in the polymorphic layer. ML: molecular layer; O-GL: outer granular layer; O-PyL: outer pyramidal layer; 1-GL: inner granular 
layer; 1-PyL: inner pyramidal layer; PoL: polymorphic layer, CC: corpus callosum. 



atrophy (reduced dendritic length) and reduced spine density (spine 
loss) (Malone et al., 2008; Martinez-Tellez et al., 2005). In their ultra- 
structural study on STZ-induced type 1 diabetes, Hernandez-Fonseca 
et al. showed that the diabetic brain is associated with degenerative 
changes of neurons and glia, presynaptic vesicle dispersion in swol- 
len axonal boutons and fragmentation of neurofilaments at 8 weeks 
(Hernandez-Fonseca et al., 2009). 

Clinical relevance of diabetes-induced cortical/striatal damage 

Frontal gyrus atrophy may contribute to cognitive dysfunction 
(Musen et al, 2006; Wessels et al., 2006, 2007), and striatopathy 
may lead to movement disorders and somatosensory abnormalities 
in clinical diabetic patients (Abe et al., 2009). It has been suggested 
that diabetic brain damage are secondary to vascular complications, 
such as increased permeability of microvasculature, high plasma os- 
molality and hypoperfusion (Jakus and Rietbrock, 2004). Interesting- 
ly, the striatum and frontal cortex, the two regions most affected by 



STZ-treatment, are also the regions known to be especially susceptible 
to ischemic damage (Lei et al., 1998; Sutherland et al., 1991 ). In addition, 
the pathohistological changes in the striatum and cortex we observed in 
the STZ-induced rats showed resemblance to what had been reported for 
rats subjected to chronic cerebral hypoperfusion (Chen et al., 2011; 
Shibata et al., 2004). On the other hand, the presence of mitochondrial 
abnormalities suggests that glucose neurotoxicity could have also played 
a role in causing diabetic brain injuries (Edwards et al., 2010; Russell 
et al., 1999). 

Limitations to the study 

There were several limitations to this study in terms of experimental 
design. 1) Ketoacidosis frequently occurs in type 1 diabetic patients 
without glycemic control, which can cause brain damage and perma- 
nent neurological deficits (Glaser et al., 2001). In this study, the level 
of ketone bodies in the plasma of the diabetic rats was not monitored. 
Possible occurrence of ketoacidosis cannot be ruled out. However, this 
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is unlikely since cerebral edema, a hallmark for ketoacidosis, was not 
observed in our diabetic rats. 2) Most of the histology and immuno- 
histochemistry techniques used in this study (i.e., silver staining, LFB 
staining and Nissl staining) are dedicated to detect demyelination 
and axonal/neuronal damage. More delicate and sensitive histology/ 
immunohistochemistry techniques may reveal more subtle and specific 
pathological changes associated with diabetic brain damage such as 
synaptic loss and abnormal microglia activity (DeVisser et al., 2011; 
Idan-Feldman et al., 2011). 3) Given the complexity of WM configura- 
tion in the striatum and cortex, and the fact that only a limited number 
of animals were subjected to histology/immunohistology evaluations, 
no attempt was made to establish semi-quantitative correlations be- 
tween the diabetes-related changes in diffusion indices and specific 
histopathological alterations. 4) Without glycemic control, the gen- 
eral conditions of the STZ-induced animals may deteriorate rapidly 
(Hoybergs et al., 2008). This is also the reason we performed DTI at 
4 weeks after induction, but not at later time points when cognitive 
impairments are shown to be prominent (Idan-Feldman et al., 201 1 ). 
Further studies are needed to establish whether or not changes of 
diffusion measures in the diabetic striatum and cortex animals are 
associated with cognitive/motor dysfunction. 5) Further studies on 
animal models of type 2 diabetes are needed, so are the studies designed 
to investigate the effect of glycemic control on the progress of brain 
complications. 

There are also some limitations with respect to the image processing 
pipeline used in this study. First, Leemans and Jones have shown that 
the b-matrix needs to be adjusted during the correction of eddy cur- 
rents and residual motion drifts (Leemans and Jones, 2009), and the in- 
tensity of the diffusion data has to be modulated with the Jacobian 
determinant due to residual phase-encoding eddy currents (Jones and 
Cercignani, 2010). Such steps were not performed in this study and 
could have potentially affected the accuracy of tensor fitting. 

Second, selection of the type (i.e., isotropic vs. anisotropic) and 
kernel size of spatial smoothing is known to influence the VBA results 
on DTI data comparisons. For example, it has been an empirical rule of 
thumb that the FWHM should be at least twice the voxel dimension 
in order to get statistically robust results (Jones et al., 2005). In addi- 
tion, anisotropic smoothing is known to provide better results than 
isotropic smoothing when performing VBA on DTI data (Van Hecke 
et al., 2010). In this study, we used isotropic smoothing with a Gaussian 
kernel having a FWHM three times the dimension of the voxel size. An- 
isotropic smoothing was not performed due to its technical complexity. 

Finally, we used affine transformation in constructing the population- 
specific template and for image spatial co-registration. The use of more 
advanced non-rigid co-registration algorithms and tensor reorientation 
strategies could potentially give better results by improving the spatial 
and orientational alignment of the diffusion tensor images (Van Hecke 
etal., 2008, 2011). 



Conclusion 

In conclusion, we showed that, without glycemic control, the STZ- 
induced type 1 diabetic rats had brain injuries observable to DTI as 
early as 4 weeks after induction. The cerebral abnormalities associated 
with STZ-induced type 1 diabetes are characterized by significantly 
reduced FA in the motor/somatosensory cortex and striatum. The 
FA changes in the diabetic cortex and striatum are likely manifesta- 
tions of demyelination and axonal degradation in these regions. 
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